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Abstract. We study the addditon problem for strongly matricially free random 
variables which generalize free random variables. Using operators of Toeplitz type, we 
derive a linearization formula for the matricial R-transform related to the associated 
convolution. It is a linear combination of Voiculescu's R-transforms in free proba- 
bility with coefficients given by internal units of the considered array of subalgebras. 
This allows us to view this formula as the matricial linearization property of the R- 
transform. Since strong matricial freeness unifies the main types of noncommutative 
independence, the matricial R-transform plays the role of a unified noncommutative 
analog of the logarithm of the Fourier transform for free, boolean, monotone, ortho- 
gonal, s-free and c-free independence. This paper treats the case of two-dimensional 
' arrays. However, all results can be generalized to arrays of any finite dimension and 

a more general version of this paper will appear elswehere. 
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1. Introduction 

In this paper we study the addition problem for strongly matricially free random 
variables [12]. These results extend Voiculescu's results [22] on the addition of free 
random variables (generalized by Maassen [15] and Bercovici and Voiculescu [2]) and 
include those for the addition of random variables associated with other fundamental 
types of noncommutative independence (monotone [16], boolean [20]), generalizations 
of these (conditionally free [4], conditionally monotone [8]) as well as those related to 
free subordination (s-free [10], orthogonal [10]). 

In classical probability, the addition problem for classically independent random vari- 
ables is related to the classical convolution of probability measures \i\ * fi 2 and to the 
^ ' logarithm of the Fourier transform which linearizes this convolution, 

% ■ 

(1.1) logF M1 ^ 2 =logF w +logF M2 , 

where is the Fourier transform of fi. However, the analogous addition problem for 
noncommutative random variables is more involved since there is no single notion of 
noncommutative independence. 

In free probability [21-25], there is a remarkable analog of the above formula, in which 
the classical convolution is replaced by the free convolution fii EE /i 2 and the role of the 
logarithm of the Fourier transform is played by the R-transform, 

(1-2) ^Wffl/U2 = ^Pl + R-w> 

where is the R-transform of \x. However, the results on additive convolutions asso- 
ciated with other types of noncommutative independence imply that the R-transform 
does not retain the linearization property in the general noncommutative framework 
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and other transforms, such as the reciprocal Cauchy transform [16], the K-transform 
[20], or the c-free R-transform [4], have to be used to describe these convolutions. 

This situation was one of our motivations to look for new types of independence which 
would unify the existing fundamental types and for which the associated transform could 
play the role of a noncommutative logarithm of the Fourier transform (some motivation 
came also from [14]). In this work, we focus our attention on a generalization of freeness 
called strong matricial freeness [12], in which arrays of non-unital *-subalgebras (Aij) 
of a given unital *-algebra A replace families of unital *-subalgebras of free probability. 
We assume, however, that each algebra Aij has an internal unit ljj which is a projection 
and that all these units, together with 1 A , generate a commutative *-subalgebra X called 
the algebra of units. Moreover, one distinguished state is replaced by an array of states 
(<Pi,j) on A, where by a state we understand a complex- valued normalized positive linear 
functional. 

We have shown in [12] that strong matricial freeness unifies the fundamental types of 
noncommutative independence. Namely, including some additional cases distinguished 
in this paper, we obtain a correspondence 

shapes of arrays ^ types of independence, 

where, in particular, square, lower-triangular and diagonal arrays correspond to free- 
ness, monotone independence and boolean independence, respectively. This scheme 
also includes the notions of s-freeness and orthogonal independence and leads to the 
correspondence 

addition of rows ^ convolutions of measures, 

where by a row we understand the corresponding sum of variables. All binary additive 
convolutions obtained in this fashion are induced by addition of strongly matricially 
free random variables. Since all convolutions considered in this paper will be additive, 
the adjective 'additive' will be usually omitted in the sequel. 

Our scheme shows that strongly matricially free random variables are basic noncom- 
mutative variables whose addition leads to well-known convolutions. Thus, it is natural 
to define the strongly matricially free convolution of the array of distributions (/^j) as 
the distribution of the sum 

(1.3) .1 V, /;? 

of strongly matricially free random variables in a distinguished state (p, where /ijj is 
the distribution of a iy j e Aij in the state ifij, denoted by EEkj/^j- If we work in the 
category of C*-algebras and the variables are self-adjoint, the considered distributions 
can be identified with compactly supported probability measures on the real line. 

For simplicity, let us specialize to the case of arrays (dij) in which (i, j) e J, where J 
is a subset of the Cartesian product {1,2} x {1,2}, often omitted in our notations. The 
symbol El used for our convolution operation does not mean that the distribution of A 
is the free convolution of an array of distributions. In particular, one cannot change 
the 'matricial order' in which the distributions appear. However, if the array (a^) is 
square and the corresponding distributions (fJ-ij) are row-identical, then 
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where /i^i = yU li2 = /ii and [i 2 ,i = ^2,2 = A*2- Therefore, in this case, the new con- 
volution can be viewed as a decomposition of the free convolution. Moreover, this 
framework also gives decompositions of other binary convolutions in terms of (j^ij) if 
we consider subarrays of square arrays or relax the assumption that the distributions 
are row-identical. In this fashion we obtain boolean, monotone, conditionally free and 
conditionally monotone convolutions, as well as s-free and orthogonal convolutions re- 
lated to the subordination property of the free convolution discovered by Voiculescu 
[25] and generalized by Biane [3]. For other results on s-free convolutions, see also [11] 
(mulitplicative case) and [17] (multivariate case). 

We would like to find a suitable noncommutative analog of the logarithm of the 
Fourier transform. For that purpose we study strongly matricially free Toeplitz opera- 
tors, similar to those in free probability studied by Voiculescu [22] and Haagerup [7] (an 
extension of which to conditional freeness has been proposed recently [18]). We follow 
the approach of Haagerup who used adjoints of the Toeplitz operators introduced in 
the original approach of Voiculescu. Namely, we take 

(1-5) <>,, • 

where £i,j,£*j are strongly matricially free creation and annihilation operators, respec- 
tively, living in the strongly matricially free Fock space and f it j is a polynomial for 
any e J. A novelty, as compared with the free case, is that the constant term of 
fi,j(£*j) is equal to the internal unit ljj multiplied by some complex number. 

This allows us to derive an equation for the Cauchy-transform Qa of some distribution 
of A in the state p with an operator-valued argument taken from the algebra of units 
X. This equation is of the same form as in the case of scalar-valued arguments, namely 

(1.6) Ga(- z +K a (z^J = z 

in some 'neighborhood' of zero. However, the argument of Qa is an invertible X-valued 
power series which has a bounded inverse, with TZa(z) playing to role of an X-valued 
analog of the R-transform. Another important point is that the distribution of A con- 
sidered here, although scalar-valued, cannot be identified with the collection of moments 
of A in the state ip, but rather with the collection of moments of A alternating with 
elements of X, as in the operator- valued free probability [23,24]. 

This distribution can be viewed as a noncommutative extension of the strongly ma- 
tricially free convolution considered above, by abuse of notation also denoted EEkj/^ij, in 
which we convolve the array (fJ-ij) in a more noncommutative fashion. Using it, we can 
prove the addition formula for X-valued R-transforms involved which not only extends 
the scalar- valued linearization formula (1.2), but also takes a nice form 

(1.7) n A = Y, n ^ 

where TZij(z) = Ri t j(z)l it j and Rij denotes the scalar- valued R-transform of //jj for 
any e J- In fact, this formula shows that convolutions associated with the main 
types of noncommutative independence can be linearized by the X-valued R-transform, 
called matricial R-transform, except that we have to use noncommutative distributions 
of A and modify the concept of linearization. Essentially, this noncommutative setting 
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turns out natural for proving existence of a noncommutative analog of (1.1) and (1.2) 
which is suitable for various types of independence. 

We can also view this formula as a linearization property of the R-transform of a new 
type, since the matricial R-transform is a linear combination of the R-transforms of 
Voiculescu with X-valued coefficients. Let us stress that this linear combination is very 
special since its coefficients are exactly the internal units of the considered subalgebras. 
We should also add that a similar formula holds for all states from the array (</>ij) 
if we take suitably truncated operators and only when considering arrays of Cauchy 
transforms of this type we can prove uniqueness of the matricial R-transform. 

In particular, if the array is square and distributions are row-identical, formula (1.7) 
gives additivity of the scalar-valued R-transforms since two different unit decomposi- 
tions lead to decompositions of the R-transforms, namely 

2 2 

for each i e {1, 2} and thus (1.2) becomes its consequence. Let us remark in this context 
that the fact that the distributions are row-identical does not mean that such are the 
states. In fact, in the GNS construction given in [12] the off-diagonal states are by 
construction different from each other and different from the diagonal states. 

Conversely, we can arrive at the matricial R-transform starting from free probability 
and using operatorial subordination [10]. For that purpose, we begin with a pair of 
free Toeplitz operators ai,a 2 e B(F(H)), where J r ('H) is the full Fock space over a 
two-dimensional Hilbert space with an orthonormal basis {ei,e2}. The results of [10] 
give decompositions 

(1.9) aj = Xj +Xj for j e {1,2} 

where Xj and Xj are equal to a,j restricted to ^(Cej) and its orthogonal complement 
J-^Cej) 1 , respectively, where j e {1, 2}. In this case, the ^-distribution of x\ agrees with 
the (/^-distribution of X 1 and the ^-distribution of x 2 agrees with the (/^-distribution 
of X 2 (these distributions form a matrix and can be denoted /ii tl , /ii i2 , £t 2 , 2 and fj, 2 ,i, 
respectively), where (pj is the state defined by the vector e,- and j e {1,2}. 

Treating the R-transforms of the distributions of a\ and a 2 as operators on J-(H) 
and decomposing the unit of S(J r ('H)) as described above, we arrive at (1.7) with 
A = ai + a 2 . Then it remains to relax the assumption that the distributions are row- 
identical to obtain the general case. Therefore, another way to arrive at the matricial 
R-transform is to follow the 'subordination path' 

( ai \ f Xl x 1 \ ( n hl n h2 \ 

\a 2 J ~* \ X 2 x 2 J ~* \ 7*2,1 ^ 2 , 2 J ' 

namely first decompose free random variables as in (1.9) and then assign to them the 
associated operatorial R-transforms induced by unit decompositions (1.8). 

If we deform these R-transforms and thus allow for different R-transforms in each 
row, the variables corresponding to rows become conditionally free with respect to 
(ip, tp), where tp is any extension of ifi and <p 2 . In this context, let us point out that the 
advantage of using strong matricial freeness instead of conditional freeness is that the 
first one includes monotone independence, whereas the second one requires additional 
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(rather restrictive) assumptions on the considered subalgebras, although satisfied in the 
case of algebras of polynomials used when studying convolutions [6]. 

At the same time, our theory has some new features since it provides a framework 
placed between scalar- and operator-valued free probability. An important property is 
that despite the fact that the internal units are not identified, the considered states 
are scalar- valued. A similar approach was used to define matricial freeness [12], which 
could also be called 'weak matricial freeness', especially when compared with strong 
matricial freeness. Addition of (weakly) matricially free random variables will be studied 
elsewhere. 

Let us remark that our construction of a noncommutative logarithm unifying the 
logarithm of the Fourier transform and the K-transform given in [9] is of a different 
nature than the approach presented in this work, although it is possible that our present 
approach can be extended in that direction. 

The remainder of this paper is organized as follows. Section 2 is devoted to the 
concept of strong matricial freeness. The corresponding convolution is studied in Section 
3. Suitable Toeplitz operators living in the strongly matricially free Fock space are 
introduced and studied in Section 4. In Section 5, we define the matricial R-transform 
and prove the corresponding linearization formula, which is the main result of this paper. 
In Section 6, we give conditions under which the matricial R-transform is unique. In 
Sections 7 and 8, we study strongly matricially free convolutions from a combinatorial 
point of view. 

2. Strong matricial freeness 

Let us recall the basic notions related to the concept of strong matricial freeness [12]. 
Let A be a unital *-algebra with an array (Aij) of non-unital *-subalgebras of A and 
let (tpij) be a family of states on A (there is some similarity to freeness with infinitely 
many states [5]). Further, we assume that each Aij has an internal unit ljj which is a 
projection for which alij = l^a = a for any a e Aij and that the unital subalgebra X 
of A generated by all internal units is commutative. 

Crucial is the fact that internal units are not identified with the unit of A and 
therefore additional conditions on moments involving these units are needed. In order 
to state these conditions, we need subsets of (/ x I) m of the form 

T m = {((«1, *2), (*2, *3), -■-;(* )) : k # 12 # • • • # i m ) 

where / is an index set and meN, with the corresponding union denoted 

00 

r = |J r m . 

m=l 

Objects (algebras and their elements, states, etc.) labelled by indices of the form 

for any j and i ¥= j for any i # j, respectively, will be called diagonal and off-diagonal. 

Important is the difference between diagonal and off-diagonal objects. 

Definition 2.1. We say that (Uj) is a strongly matricially free array of units associated 
with (Aij) and (<fij) if for any diagonal state ip it holds that 

(a) tp(u\au2) = ip(ui)ip(a)ip(u2) for any a e A and u±, it 2 e X, 

(b) (p(l id ) = 5i d for any i,j, 
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(c) if a k e A ik>jk n Ker<^, 



,]ki 



where 1 < k ^ m, then 

cp(aa 2 ...a n ) if ((i 1 J 1 ),...,(i m ,j m ))eT 
otherwise 



ip(al iujl a 2 ...dm) = 



where a e A is arbitrary and («i, ji) (i m ,jm)- 

The main condition of strong matricial freeness reminds freeness as the definition 
gives below shows. 

Definition 2.2. We say that *-subalgebras (Aij) are strongly matricially free with 
respect to (<Pij) if the array of internal units (hj) is the associated strongly matricially 
free array of units and 

(2.1) ^( a i«2 • • • a n ) = whenever a k e A ikJk n Kery? ife)jfc 

for any diagonal state ip, where ¥=■■■¥= (i n ,jn)- 

Using these conditions, one can easily show that multiplication of weak and strong 
matricially free random variables in the kernel form reminds multiplication of matrices 
and arrays of units remind matrix units, except that in the strong case the elements 
which are on their diagonals survive only in the last matrix. Naturally, the array of 
variables (aij) in a unital *-algebra A is called strongly matricially free with respect to 
(<Pij) if there exists an array of projections (lij) which is a strongly matricially free array 
of units associated with A and (ifij) and such that the array (Aij) of *-subalgebras of 
the form Ai j = alg(aij, a*j, lij) is strongly matricially free with respect to ((fiij). If A 
is a C*-algebra, we assume that the subalgebras Aij and X are C*-subalgebras of A. 

We shall consider an array of states (<Pij) on A in which diagonal states coincide with 
a distinguished state ip, whereas the off-diagonal ones are given by tpij = tpj, where 
i j and 



states (they were called 'conditions' in [12]). In particular, this implies the normalization 
conditions ^-(1^) = for any i,j, k. 

Finally, if the distribution of a^j in the state ipij does not depend on j, we will say 
that the array (ojj) is row-identically distributed. 

Proposition 2.1. Let (Aij) be strongly matricially free with respect to (<fi,j), where 
tfjj = (p for any j and tpij = tpj for any i ¥= j ■ Let us consider the mixed moment 
Vi,j(aia 2 . . . a n ), where a k e A ik:jk and (hj^ ¥=...¥= (i n ,jn)- 

(1) If a n is diagonal, j n # j and the state ipij is off-diagonal, then the moment 
vanishes. In particular, the off-diagonal state (p it j vanishes on any A k ,k for 
k^j. 

(2) // a n is off-diagonal and the state Lp i; j is diagonal, then the moment vanishes. 
In particular, a diagonal state vanishes on any off-diagonal subalgebra. 

(3) If a k is diagonal for some 1 ^ k <n and a r e KeTtp ir j r for k < r ^ n, then the 
moment vanishes. 




Proof. These are straightforward consequences of the definition of strong matricial 
freeness. ■ 
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Definition 2.3. By the strongly matricially free Fock space over the array (Hij) we 
understand the Hilbert space of the form 

oo 

(2.3) A/- = Cfi0 © KZ®Kl2®---®K?:,r m , 

m=l i 1 ^...^ti m 
n 1 ,...,n m GN 

where Q is a unit vector, with the canonical inner product. 

In particular, when the array is square and consists of one-dimensional Hilbert 

spaces Hi j = Ce^j for i,j e {1,2}, we have 

oo 

M = ©A/" (m) , 
where the first few summands are of the form 

a/" (0) = cn 

AA 1 ) = C eiil 0Ce 2)2 

A/" (2) = Ceg©Ceg©C(e li2 ®e 2i2 )eC(e 2il ® eiil ) 

A/" (3) = Ce®?©Ceg©C(e 2il (g)e®?)©C(e 1)2 (g)eg)©C(eg(g)e lil ) 

© C(eg ® e 2 , 2 ) © C(ei )2 ® e 2) i ® ei,i) © C(e 2 ,i ® ei )2 ® e 2)2 ), 

etc. Note that, in contrast to the (weakly) matricially free Fock space, diagonal vectors 
appear only at the ends of the tensor products. 

Let us distinguish the closed subspaces Af it j spanned by simple tensors which begin 
with e it j, where e J. Then we have a decomposition 

(2.4) A/" = Cfi©0A/;,, 

and the diagonal unit ljj is the projection onto CQ ©A/}j = J-(Cej t j), whereas the 
off-diagonal unit ljj is the projection onto A/"© J-^Ce^j). 

Remark 2.1. Let p be the projection onto Cf2 and let pij be a projection onto Afij for 
any i,j. Note that we have 'orthogonal decompositions' of the unit from A = B(J\f), 

(2.5) 1 A = li,i + li, 2 = 1 2 ,2 + l 2 ,i, 
as well as 

( 2 - 6 ) U =P + 2p« = 2*J' 

where g^i = p, gi j2 = p 2i2 , q 2i i = and g 2j2 = pi ;2 + p 2 ,i, which will be of use in our 
study of the matricial R-transform. 

Definition 2.4. Let A = (aiij) be an array of positive real numbers and let (Hi,j) = 
(Ceij) be the associated array of Hilbert spaces. By the strongly matricially free creation 
operators associated with A we understand operators of the form 

(2.7) £ id = a h] a*t{e hj )o-, 

where a : Af — > ^*(©j is the canonical embedding in the given free Fock space 
and the £(ejj)'s denote the canonical free creation operators. By the strongly matricially 
free annihilation operators we understand their adjoints. 
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Proposition 2.2. The array (Aij) of *-subalgebras of B(M), where A id = alg(i^j^i*j) 
for any e J, is strongly matricially free with respect to (fij)- 

Proof. The proof is left to the reader since it is very similar to that in the free case 
(see also [13, Proposition 4.1]). Let us only observe that the relation 

(2-8) = oZjlu 

implies that e Aij since a i: j is assumed to be positive. ■ 

Remark 2.2. In a similar way one defines the concept of matricial freeness, which, 
especially when compared with the strong case, can also be called 'weak matricial 
freeness'. The difference is that in the definition of the array of matricially free array 
of units we replace the sets T n by slightly larger sets 

^2) J2 ) j ■ ■ ■ j [}mi jm)) ■ (*1,*2) (i2,is) ^ ■■■ ^ (i mi im+l)} j 

where all 'matricially related' pairs of indices are allowed. In this case, in order to 
obtain a non-zero mixed moment of type v 9 (°i a 2 • • • Q> n ) it is also necessary that a n is 
diagonal, but diagonal variables can also appear in the middle of non-zero moments of 
the type considered in Proposition 2.1. 

3. Convolutions 

From the results of [12] it follows that different shapes of (a it j) correspond to different 
types of noncommutative independence. The convolution which unifies the associated 
convolutions is the convolution of an array of distributions defined below. 

Definition 3.1. Let (a^j) be a two-dimensional array of strongly matricially free ran- 
dom variables with the corresponding array of distributions (faj) in the states (fij)- 
The distribution of the sum 



(3.1) A = J> 



1,3 

i,3 



in the state <p is called the strongly matricially free convolution and will be denoted 
£Bi,jfJ>i,j by analogy with the free convolution. In these notations it is understood that 
the summation runs over pairs (i,j) e J, which is often omitted. 

In particular, the strongly matricially free convolution generalizes the free convolu- 
tion, which justifies our notation. For simplicity, let us consider a two-dimensional array 
(a it j) of random variables and an array of states given by 



(3.2) 



<Pl,l <Pl,2 \ = I <P 
¥2,1 ^2,2 J \¥l ¥ 

where ip is a distinguished state on A and v ? i)V 9 2 are conjugate states, with respect to 
which the considered array is strongly matricially free. 

For the sake of greater generality, by a two-dimensional array we shall understand a 
subbarray of a square two-dimensional array. Consider sums 

(3.3) Oi = a id 

j 

where % = 1,2 and it is understood that the summation runs over those j's, for which 
(i,j) e J. In other words, we add variables which appear in each row separately. 
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Clearly, these sums may reduce to one term if only one variable appears in a given 
row. Let us remark that the array is not assumed to contain the diagonal as in the 
case of limit theorems [12,13]. This allows us to consider arrays corresponding to s-free 
independence and orthogonal independence. 

In the theorem given below we state explicitly how addition of rows corresponds to 
convolutions associated with various notions of noncommutative independence and how 
these convolutions can be decomposed using the strongly matricially free convolution. 
By the (^-distributions of A and ai we understand the collections of moments of A and 
Oj, respectively, in the state ip. Each pair of variables is denoted {ai, a 2 }, but one should 
remember that in some cases (3-5) order in which they appear is relevant. 

Theorem 3.1. Assume that (a^j) is strongly matricially free with respect to ((pij) and 
row-identically distributed with distributions (faj). Let fj, and fa be the tp- distributions 
of A and a i7 where i e {1,2}, respectively. 

(1) If the array is square, then {a\, a 2 } is free with respect to ip and fi = fa EB fa ■ 

(2) If the array is diagonal, then {oi,a 2 } is boolean independent with respect to p 
and [i = fa ty fa. 

(3) If the array is lower-triangular, then {ai,a 2 } is monotone independent with re- 
spect to p and /i = fa o fa. 

(4) If the array is upper-anti-triangular, then {ai,a 2 } is s-free independent with 
respect to (p, p\) and /i = fa □ fa. 

(5) If the array consist of one column, then {fa, a 2 } is orthogonally independent with 
respect to (p, p\) and fi = fa \- fa- 
Proof. Since A = fa + a 2 , cases (l)-(3) are immediate consequences of [12, Propo- 
sition 4.1], where it was shown that {ai,a 2 } is free, boolean independent or monotone 
independent, depending on whether the array is square, diagonal or lower-triangular, 
respectively. In order to show cases (4) and (5) it suffices to show that {ai 5 i + ai )2 , «2,i} 
is s-free and {ai 5 i,a 2j i} is orthogonal under {ip, pi). However, these properties follow 
from the GNS representation in the subordination context given in [10]. ■ 

We have distinguished only those convolutions which are non-trivial and cannot be 
reduced to the other ones. This eliminates upper-triangular arrays (formally they cor- 
respond to anti- monotone independence and its conditional counterpart), which can be 
reduced to lower-triangular ones, as well as lower- anti-triangular arrays (i.e. those in 
which there is no a^i), which can be reduced to upper- anti-triangular ones (i.e those in 
which there is no a 2)2 ). Finally, we have omitted arrays consisting of one row, in which 
case the (/^-distribution of A agrees with the (^-distribution of the diagonal variable a^i 
and the (/^-distribution of A agrees with that of the off-diagonal variable ai i2 (all mixed 
moments of two variables vanish). 

If we lift the assumption that the variables are row-identically distributed, we obtain 
the conditionally free additive convolution. In a similar manner, a generalization of the 
case of lower-triangular arrays leads to conditionally monotone variables of Hasebe [8]. 

Theorem 3.2. If is a square array of random variables in A which is strongly 
matricially free with respect to (pij) and has distributions (faj), then 



(3.4) 



fflij fkj = Vl,2) ffl (f*2,2, A*2,l) 
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and {ai,a 2 } is conditionally free with respect to (<p,ip), where ip is any state on A 
extending </?iU 2 an d V^Ui, where Ai = alg(Ai,i, Ai, 2 ) and A2 = alg(A 2> i, A 2 , 2 ). 

Proof. It suffices to prove that variables of the form 

a° = a* - ip(ai) = J]( a i,j ~ V'Kj)) 
j 

satisfy the equation of c-freeness with respect to tp, namely 

^(<<---<)=^(<M<)--^(aD 

where i\ i 2 ¥=■■■¥= i n and ip is any extension of <£>i|.a 2 and <p 2 \ai- The assumption 
on if; gives 

since y?(ai )2 ) = ^(^2,1) = as well as ^2(01,1) = V?i( a 2,2) = by Proposition 2.1. 
Therefore, our variables can be decomposed as 



a 







a' + a" + 



where the first two terms belong to diagonal subalgebras, namely 

a- = a M - ^(a^jly 

for each z and 

a" = (</?(°i,i) - ^(ai )2 ))li,i, a" = (<p(a 2l2 ) ~ ip(a 2 ,i))l 2 , 2 
whereas the third one is an element of an off-diagonal one, i.e. 

61 = 01,2 - ip(a lj2 )l li2 and b 2 = a 2A - ^(a 2 ,i)l 2 ,i- 
Now, using Proposition 2.1(2) again, we obtain 

vK< ■■■<) = <p«< ■ ■ ■ <) + vK< ■■■<)• 

Next, we claim that 

¥»(«...<,) =0. 

Namely, in view of Proposition 2.1(3), the diagonal variables can appear only as the 
last ones in non-zero mixed moments of kernel form. This implies, however, that 

<fi(a° h a% ...a' in ) = tp(b h b i2 . . . b^XJ = 

by strong matricial freeness (note that this kind of reduction is not possible in the case 
of weak matricial freeness). It remains to use property (a) of Definition 2.1 to write 

v(« ■■■<) = - ^wJM<< • • • O 

and use an inductive argument to finish the proof. ■ 
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4. TOEPLITZ OPERATORS 

In order to study the addition of free random variables, Voiculescu used Toeplitz 
operators [22]. They turn out to be closely related to R-transforms and can be used 
to prove that they are additive under free convolutions. A new proof of additivity was 
presented by Haagerup [7] who used the adjoints of Toeplitz operators of the form 

(4.1) a = £ 1 + f(t 1 ) and b = £ 2 + g(£*) 

where £i,£ 2 are isometries on the full Fock space J-'('H), where "H is a two-dimensional 
Hilbert space with orthonormal basis {ei,e 2 }, given by tensoring on the left by e\ and 
e 2 , respectively, and where /, g are polynomials. Using only Fock space techniques, he 
derived the addition formula for the R-transforms. 

We would like to generalize this approach to the case of strong matricial freeness in 
order to find a suitable analog of the R-transform. 

Definition 4.1. Let (£ij) be the array of strongly matricially free creation operators 
on M and let fij be a polynomial for any e J. Operators of the form 

(4-2) ",, ■ LA'*,) 

where e J and the constant term of fij is the internal unit multiplied by a 
complex number, will be called strongly matricially free Toeplitz operators. 

These operators are similar to free Toeplitz operators used by Voiculescu and Haagerup. 
We will need the distributions of these operators in the states (tpij) defined by the array 
of unit vectors (i\j), where 

Qjj = Q and Qij = ejj for % # j, 

which replace the single vacuum vector which suffices in the free case. It is easy to see 
that the restriction of a i: j to the closed subspace of M spanned by vectors 

has the Toeplitz form for any e J. Moreover, in the case of row-identically 

distributed square arrays these operators are closely related to those in free probability 
since the sum of those lying in the j'-th row give a decomposition of the j-th free Toeplitz 
operator (see Introduction). 

Proposition 4.1. The R-transform of the distribution fiij of the operator in the 
state <fij is given by 

(4-3) H, ..,(:) f,M,:). 

where e J and the constant term of f\j is a complex number. 

Proof. The proof is similar to that in [7], but some modifications need to be intro- 
duced. For fixed e J, take the vector 

oo 

/',,(-- ) = (1 " zki)- 1 ^ = fly + 2 z n al 3 ef- ® D,, 

n—l 

where we set e®" (x) Q = e®" for any j and assume that \(Xijz\ < 1 in order to get a 
convergent series. Then 

kjPi,j( z ) = ~(Pi,j( z ) -Vij) 
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for < \z\ < \otij\ 1 . Moreover, since £*j^i,j = an d 

^W® n «)=«i J -(eg B - 1) ®nu), 

we can see that 

?i,jPi,j( z ) = za i,jPi,j( z ) 

and thus 

<*jPiA z ) = -(piA z ) - + fi,j( a lj z )Pi,j( z ), 

which leads to the equation 

+ fij( a h z ) - a ijj = 

Now, if |z| is sufficiently small and positive, the operator on the left hand side is 
invertible, which gives 

Gij + fijialjz^j = z, 

where Gij denotes the Cauchy transform of the (^-distribution of a,ij, and that im- 
plies that its R-transform has the desired form by the uniqueness of the R-transform. ■ 

It should be noted that in the study of (^-distributions of the strongly matricially free 
convolution it suffices to use strongly matricially free Toeplitz operators (a^) and their 
sum A since, in view of Proposition 2.2, the array (ajj) is strongly matricially free with 
respect to (<fij) and this implies [12, Propositions 2.2-2.3] that the mixed moments of 
these variables are uniquely determined by the marginal distributions and thus by the 
corresponding R-transforms (Rij)- Therefore, if we start with an arbitrary array of 
strongly matricially free random variables and the same distributions as those of (cbij), 
then we obtain the same mixed moments. Thus, without loss of generality, we can use 
strongly matricially free Toeplitz operators with given R-transforms. 

It can also be justifed that in this study one can use polynomials as R-transforms 
instead of infinite series 

oo 
n—l 

for any e J- This is because any mixed moment of order m of the sum of arbitrary 
strongly matricially free random variables can be expressed in terms of a finite number 
of coefficients ^(n) for 1 ^ n ^ m. 

These results give us a hint that one might be able to use strongly matricially free 
Toeplitz operators to obtain an analog of the R-transform. Following Haagerup's ap- 
proach [7], we shall first define suitable vectors of geometric series type and then act 
on these vectors with Toeplitz operators. 

Let (£ij) be the array of strongly matricially free creation operators living in the 
strongly matricially free Fock space M and let 

be the corresponding sum of creation operators. Consider the vector 

p(z) = (l-zL)- 1 ^, 
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where z is a complex number such that \z\ < Q]ij l^.il 2 ) -1 ; which ensures convergence 
of the associated series. 

Lemma 4.1. The sum A = Yjij a i,j °f strongly matricially free Toeplitz operators 
satisfies the equation 

(4-4) Ap(z) = -(p(z) - n) + J f^al^hjpiz) 

where p(z) = (1 — zL)~ l Q and < \z\ < ■ la^l 2 ) -1 . 
Proof. Clearly, as in the free case, we have 

L P (z) = -( P (z)-n). 

z 

Let us now decompose the vector p(z) into sums of two vectors in two different ways. 
These vectors will be eigenvectors of annihilation operators from each row of the array 
(£fj). Thus, using the same notation as in the proof of Proposition 4.1, we get the 
action of the diagonal annihilation operators 

e hPjA z ) = a h z PjA z ) 

for j g {1,2}. Moreover, for the off-diagonal ones we have 
where i j and 

771,2(2) = P(z)-Pi,i(z) 

V2,l(z) = p(z) - p 2 , 2 (z). 

Let us show the first of these equations. We have 

00 

r 1>2 ( P (z) - Pl>1 (z)) = 2^^iV B -^,i)n 

71=1 

00 

= 2^/1,2^"^ 

71=1 

00 



71=1 

,2 



= a li2 z(p(z) -pi,i(z)) 

since ^* 2^1,2 = ot\ 2 li,2 and li^-k"' -1 ^ = L n_1 r2 — ef^ _1 \ which proves the desired 
equation. Finally, 

*lVl,l(*) = 4*,lP2,2^) = ^,1^1,2 (*) = ^,2^,1 (*) = 0. 

Therefore, 

(/2,2^2,2) +/2,l(^2,l))PW = /2,2 (tt2, 2 ^)p2,2(20 + f2,l{a\l Z )V2,l{z), 

which gives (4.4) and thus completes the proof. ■ 
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Lemma 4.2. Let tp be the state associated with the vacuum vector fL Then there exists 
e such that 

(4.5) * = V^+%Ru(z)lij-Aj j 

whenever < \z\ < e. 

Proof. Rewriting the result of Lemma 4.1, we obtain 

For small and positive \z\, the operator on the LHS is invertible and then 

WW = + Yj fiA a h z )hj -Aj to 

which, with the use of (zp(z), f2) = z and Proposition 4.1, finishes the proof. ■ 

The above lemma indicates that a natural analog of the R-transform could have a 
more noncommutative character. This is really the case and, as we argue in the next 
section, in order to obtain the corresponding linearization formula, we need to use more 
noncommutative distributions of A than those studied in Section 3. The latter are then 
obtained by restriction to the algebra of polynomials in A. 

5. Matricial R-transform 

In view of Lemma 4.2 and its similarity to the scalar- valued counterpart involving 
the scalar-valued R-transform, we are ready to define an object called an 'operatorial 
R-transform'. For that purpose we shall use the notion of a distribution of A in the 
state tp similar to that in the operator- valued free probability [23]. 

Let a be an element of a unital complex Banach algebra A with a unital closed 
subalgebra X and let tp be a normalized linear functional on A. The algebra freely 
generated by X and an indeterminate X will be denoted by X(X). The linear functional 
fx : X(X) — > C defined by fx = tp o r, where r : X(X) — > A is the unique homomorphism 
such that r(b) = b for b e X and t(X) = a will be called the distribution of a e A in 
the state tp. Quantities of the form 

tp(b ni ab n2 . . . &n m _i^n m ), where b Hk e X for 1 s$ k s$ m and meN 

will be called the Z-moments of o. In particular, the collection of moments 

&(b(ab) n ) : n e N u {0}} 

will be called the b- distribution of a. Note that we do not assume, as in the operator- 
valued free probability [23], where tp is a conditional expectation, that tp(b) = b for any 
h c; X 

It is clear that if b e X is invertible and || 6 _1 ||<|| a then the inverse of b — a 
exists and takes the form of a series 

oo 

(b-a)- 1 = Y i b~ 1 (ab- 1 ) n 

n=0 
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which converges in the norm topology. This leads to the operatorial analog of the 
Cauchy transform of the form 

00 

(5.1) g a (b) = 2 l <p{b- 1 (ab- 1 ) n ), 

due to continuity of <p, which plays the role of the Cauchy transform of the 6-distribution 
of A in the state (p. If A is a C*-algebra, continuity follows from the fact that ip is 
a state. Clearly, in this case b, b~ l and a become bounded operators on some Hilbert 
space H. 

Definition 5.1. Under the above assumptions, let /i denote the distribution of a e A 
in the state ip. If there exists an X-valued power series of the form 



n-l 

i.6 i 



(5.2) n a (z) = YiCr, 

71=1 

where c n e X for all n e N and z e C, which is convergent in the norm topology for 
sufficiently small \z\, and for which it holds that 

(5.3) Q a (\+n a {z) 

whenever \z\ is sufficiently small and positive, it will be called an operatorial R-transform 
of the distribution /j,. 

Remark 5.1. Let A be the C*-subalgebra generated by the Toeplitz operator ay and 
the unit ly and let X = C[ly], where the pair (i,j) e J is fixed. Then 

for small and positive \z\, where !Z aij (z) = Rij(z)lij and Rij is given by Proposition 
4.1, thus TZ aiJ becomes an operatorial R-transform of the </?jj-distribution of ay. 

On the level of formal power (Laurent) series, invertibility of the argument of the 
Cauchy transform of the form (5.3) is shown as in the case of scalar-valued arguments. 
Moreover, if zIZa(z) is a contraction, the inverse converges in the norm topology. These 
facts are elementary and we present them without a proof. 

Proposition 5.1. Formal power series C(z) = 2n=o Cm2;n 1 an< ^ ^( z ) = H^o^ z " +1 » 
where c n , b n e X for any n and Co = bo = 1, are multiplicative inverses if and only if 

m 

b m = Y i (- l Y 2 c «i--- c «p 

p=l n\ + ...+n p =m 

m 

Cm = 2(-i) p 2 K-K 

p=l ni+...+n p =m 

where m, ni, . . . ,n p e N. If || zC(z) — 1 ||< 1, then the power series B{z) converges 
in the norm topology. If \\ B(z) — z ||< 1, then the series C(z) converges in the norm 
topology. 
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Note that if lZ a is an operatorial R-transform associated with Q a , then the operatorial 
argument of Q a is of the form 

(5.4) Ca {z) = - z +n a (z) 

and plays the role analogous to that of the right composition inverse of the Cauchy 
transform. However, since 7Z a is operator-valued, it is not unique. 

Definition 5.2. Let A e Abe the sum of random variables (a^) in a unital complex 
C*-algebra A which are strongly matricially free with respect to (<£>»,_,-) and let X be 
its unital C*-subalgebra generated by the internal units. If an X-valued operatorial 
R-transform TZa of the 92-distribution of A takes the form 

(5.5) n A (z) = J]n i , j (z) 

where IZij(z) = Ri t j(z)l it j and Rij is the R-transform of //jj, the distribution of a,ij in 
the state <pij, it will be called a matricial R-transform of the ^-distribution of A. 

If TZa is a matricial R-transform, then the multiplicative inverse of C\ can be ex- 
pressed in terms of R-transforms of free convolutions of (faj), with orthogonal projec- 
tions (qij) introduced in Remark 2.1 being the coefficients. In particular, this is the 
case when (a i; j) is an array of strongly matricially free Toeplitz operators, but the result 
is more general. 

Proposition 5.2. If (ctjj) is the array of Toeplitz operators (4-2) and TZa is the cor- 
responding matricial R-transform, then the multiplicative inverse of Ca{z) takes the 
form 

= S (rri^j) 

where \z\ is sufficiently small and positive and (Qij) is the array of R-transforms of 
free convolutions of the distributions of (aij) in the states ((fij). 

Proof. Decomposing the internal units in terms of (qij) defined in Remark 2.1, 
we arrive at the orthogonal decomposition 

where (Qij) is the array of the form 

Qi,i Qi,2 \ _ 1 Ri,i + R22 Ri,i + R2,i 

Q2,l Q2,2 J \ R2,2 + Rl,2 Rl,2 + -^2,1 

and thus, by additivity of the R-transform, it is the array of R-transforms of free 
convolutions of ), the distributions of (ojj) in the states (tpi,j), respectively, namely 

A*l,l EE A*2,2 ^1,1 EH A*2,l 

/•*2,2 EH A*l,2 A*l,2 EH A*2,l 

This allows us to find the explicit form of the inverse of Ca(z) e B(J\f), namely 
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which converges in the norm topology for sufficiently small and positive \z\ since each 
Qij is analytic in some neighborhood of zero, as claimed. ■ 

Remark 5.2. It is not hard to show that Proposition 5.2 can be generalized to the 
case when (ajj) is an arbitrary array of strongly matricially free variables from a unital 
complex C*-algebra A. In that case, the projections (qij) are to be understood as 
suitable elements of X expressed in terms of (Uj) by the same formulas as in the case 
of Af. 

Lemma 5.1. An X-valued power series TZ(z) = Cn^" -1 converging in the norm 

topology in a neighborhood of zero is an operatorial R-transform of the (p- distribution 
of A if and only if 

m 

(5.6) J] X ip(b ni Ab n2 ...b nk _ 1 Ab nk )=0 

k=l n\ + ...+n^=m—k 

for all m^z2, where we assume that ni, . . . , nk are non-negative integers and where the 
series B(z) = 2n=o b n z n+l is the multiplicative inverse of C(z) = 1/z + lZ(z). 

Proof. In (5.1) we substitute the explicit form for the inverse of b = C(z) given 
by the series B(z) of Proposition 5.1. Note that if \z\ is small, then the condition 
|| zC{z) — 1 || = || zlZ(z) || < 1 is satisfied. Moreover, if \z\ is sufficiently small, the norm 
of B(z) can be made smaller than || A which gives absolute convergence of the 
series (5.1) in some neighborhood of zero. By continuity of ip this leads to the identity 

oo oo 



v (b ni )z n > +1 + J] v(b ni Ab n2 ) 

ni=0 rti,rt2=0 

oo 

+ J] V {b ni Ab n2 Ab n3 )z n ^ +n ^ + . . . 



ni,n2,n3=0 



= EX X ^{b nx Ab n2 ...b nk _,Ab nk )z™ 

m=l k=l n-L + ...+rik=m—k 

for small \z\, which gives necessary and sufficient conditions for 1Z(z) to be an operatorial 
R-transform of the ^-distribution of A, namely 



m 



^ J] v(b ni Ab n2 . . . b nh _ x Ab nk ) = 

fc=l ni +...+rik=m— k 

for all natural m, where we assume that the summation indices ni,...,n p are non- 
negative integers, which completes the proof. ■ 

Remark 5.3. It is easy to see that (5.6) is a recursion since it is equivalent to 

m 

<p(bm-i) = - J] J] ip(b ni Ab n2 . . . b nk _ x Ab nk ) 

k=2 ni + • • .+ri) e =m—k 

for m ^ 2 and bo = 1, and it is clear that it always has a solution. The conditions for 
the first few coefficients b m take the form 

<p(h) = -<p(A) 
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cp(b 2 ) = -ifiAh) - cpihA) - cp(A 2 ) 

ip(b 3 ) = -ip(Ab 2 ) - ifihAh) - <p(b 2 A) - </?(A 2 fei) - ipihA 2 ) - ip(AhA) - <p(A 3 ) 

etc., with bo = 1. In the scalar- valued case, we obtain a unique solution for all coeffi- 
cients b m which are expressed in terms of (^-moments of A. In the operator-valued case, 
conjugate states will be used to address the uniqueness problem in Section 6. 

We are ready to state and prove the main result of this paper, the addition formula 
for the matricial R-transform, which can also be viewed as the matricial linearization 
property of the R-transform. 

Theorem 5.1. // (aij) is an array of random variables from a unital complex C*- 
algebra A which is strongly matricially free with respect to (ipij) and (Rij) is the cor- 
responding array of R-transforms , then 

(5.7) K A (z)=J]K id (z), 

where A = 2jj a j,i an d ^-i,j{ z ) = ^i,j( z )^-i,j f or an V (hj) G J> w ^ sufficiently small 
\z\, is an operatorial R-transform of the if- distribution of A. 

Proof. If (a i; j) is the array of Toeplitz operators, the assertion is a consequence of 
Lemma 4.2. It remains to be shown that the case of arbitrary arrays (a^j) reduces to 
that of Toeplitz operators. Suppose (a^ ) is an array of elements of a C*-algebra A which 
is strongly matricially free under (ipij) with (ljj) being the associated array of internal 
units and let (Rij) be the array of R-transforms of the corresponding distributions. Let 
Rij be analytic for \z\ < e it j, where (i,j) e J and each e i; j is a positive number. Then 
the series 

K A (z)=2 l R iJ (z)l i j = J Vr,,( ( i)l, z^- 1 

hO n = l \i,3 J 

converges in the norm topology to an element of X for \z\ < e, where e = minjj e^j. For 
each natural p, consider its truncation 

K A(P) (z)=Y,K ( $(z), 

where 



= R%(z)h 3 and R^(z) = £ r,»:" 

n=l 



for any (i, j) e J and \z\ < e, and where A^> = ^ . a^j is the sum of strongly matricially 
free Toeplitz operators with .R-transforms (R\ P j), respectively. By Lemma 4.2, we have 

Qa(p) { 1" K A ( p )(z) 

\ Z 

for any p and < \z\ < 5 ^ e and some 5. Therefore, by Lemma 5.1, 



m 



fc=l n\ + ...+nk=m— k 



2 fiJ^A^b^ ■ ■ ■ bi^AtobV) = 
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for all m ^ 2, with Oq = 1. Now, by Proposition 5.2, J 's are the coefficients of the 
power series 




where Qf - is the truncation of Q it j to the polynomial of order p — 1 and is the 
array of projections of Remark 2.1. On the other hand, 6 n 's are the coefficients of the 
power series 

BW -2(rri^w)*" 

where, by abuse of notation, we use the same symbols for which are now the 

projections from X (expressed in terms of (Uj) by suitable equations, the same as in 
the case of AT). It can be seen from Proposition 5.1 or directly from the above formulas 

that bffl = b n whenever ^ n ^ p. Moreover, similar arguments as those concerning 
the moments of A in the state tp can be used for the alternating mixed moments of A 
(since 6 n 's and are linear combinations of internal units) to conclude that 

rtbVAtobjg . . . b^_AW£) = V {b ni Ab n2 . . . b^Ab^) 

whenever n\ + . . . + Uk + k — 1 ^ p. In this fashion, condition (5.6) for any given m 
can be shown to be satisfied by taking sufficiently large p. This implies that the se- 
ries TZa is an operatorial R-transform associated with Qa, which completes the proof. ■ 

In contrast to highly non-unique operatorial R-transforms of Definition 5.1, the ma- 
tricial R-transform can be made unique if we consider the whole array of states (<fi,j) 
and not only the distinguished state ip. Before we discuss the uniqueness problem, we 
shall make a few simple remarks related to the cases discussed in Theorems 3.1-3.2. 

(1) If the array is square and row-identically distributed, then the matricial R- 
transform associated with Qa takes the form 

TZa(z) = R^(z)l A + R^(z)l A 

due to unit decompositions (1.8) and can be identfied with the scalar- valued 
R-transform of //i E3 A*2- 

(2) If the array is diagonal, then the matricial R-transform associated with Qa takes 
the form 

K A (Z) = ^(^11,1+^(^)12,2 

which linearizes the extended boolean convolution. 

(3) If the array is lower-triangular and row-identically distributed, then the matricial 
R-transform associated with Qa takes the form 

TZa(z) =^(^11,1+^(^)1^ 

which linearizes the extended monotone convolution. 

(4) In the general case and with the use of notations used for c-free convolutions, 
we can write the matricial R-transform associated with Q A as 

n A (z) = Y i QiA^ 
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where Qij(z) = Qij(z)qij for any i,j and (Qij) is the array of R-transforms of 
free convolutions 

/ /ii EB Hi /Ui EB V2 \ 

V H2 S V\ V X EB ZV 2 J ' 

obtained by rewriting the array given in the proof of Proposition 5.2. 

Without any further assumptions, the matricial R-transform in the state ip is not 
unique. In most cases this fact is easy to see since the R-transform of A, say Ra, 
always exists and the corresponding scalar- valued B(z) satisfies conditions (5.6), but 
it is usually different than TZa given by (5.7). However, if the array is square and 
distributions are row- identical, then TZa agrees with Ra since = f° r an Y 

i e {1,2} (then, to show non-uniqueness, one needs to find another example). 



6. Uniqueness 

Let us investigate the uniqueness problem for the operatorial R-transform. For that 
purpose we shall study higher order states, including the conjugate states. 

Anticipating that some essential work will have to be done for operators of Toeplitz 
type, we first introduce vectors in the Fock space M which are similar to u(z), but 
which are generated by vectors e^i and e 2j 2- Let 

(6.1) Uj (z) = (l-zL)- 1 e jJ 

for any j e {1,2}, where \z\ is sufficiently small to ensure convergence. 

Proposition 6.1. Let f be a polynomial and let the constant term of f((*j) be a con- 
stant multiplied by ljj for any (i,j) e J. Then 

(6.2) /(*;>,-(*) = f(«jjz)l jtjUjd (z) + ajJ Df(a<jjz)n 

for any j e {1,2}, where Df(z) = Af(z)/z, with A/(z) = f(z) — /(0). Moreover, 

( 6 - 3 ) /(^jHW = fiA a h z ) l i^k{z) 

whenever i j and k is arbitrary. 

Proof. The proof is similar to that of Lemma 4.1. For instance, 

tlM*) = q j e, J + z£y jJ {l + zL + Z ^ + ..^ 
= a jd n + zoijjljju^z) 
for any j e {1,2}, which leads to the recursion 

{t h3 fuj 3 {z) = o^n^y V. + (zal j ri. j u j (z) 
for any natural n since £*j commutes with ljj. This, in turn, gives 

/(*;>,-(*) = f(^)h^(z) + n,,/;/(nj. ; ;)0 

for any j and any polynomial /, which proves (6.2). The proof of (6.3) is exactly the 
same as that of Lemma 4.1. ■ 
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In addition to the Cauchy transform studied in Section 5, we shall now study the 
Cauchy transforms of noncommutative distributions of a e A in the conjugate states 
ipi and y?2- This leads to the array 

oo 

(6.4) Qp) (6) = ^.((ft _ a)- 1 ) = £ ^(^(or 1 )") 

n=0 

whenever (i,j) e J and 6 e X is invertible with || b~ l ||<|| a In particular, we 

would like to study the case when a is of the form 

(6.5) Aij := PijAPij, where e J 

and A = Yjij a>i,j, with (Pij) being the array of projections given by 

(6.6) P { U-Ji,il2,2 jf » = j 

Since these projections belong to X, the operators Aij belong to A for any e J. 
A more intuitive understanding of these notions can be acquired in the Fock-space 
context, to which we return below. 

Lemma 6.1. Let A be the sum of strongly matricially free Toeplitz operators in Af with 
the matricial R-transform 1Z A given by (5.7). Then 

(6.10) Qij ^ — ^1Za(z)^ = z whenever (i,j) e J 

for sufficiently small \z\ > 0, where Qij is a short-hand notation for 

Proof. Using Proposition 6.1, we obtain 

AuJij(z) = -(uJij(z) - Ci,j) + n A (z)u itj (z) + n, ; /;/, ; (o7 ; ;)0 ; ; 

for any 6 J. This equation is very similar to (4.4), except for the term involving 
the difference quotient. However, the projection P it j maps the vector Qjj to zero and 
thus 

Aiju i:j (z) = -(uij(z) - +n A (z)u iJ (z) 
This leads to the equation 

zuij(z) = + K A (z) - Aij^J Qj 

for small \z\, as in the proof of Lemma 4.2, which proves (6.10) since (ooij(z), Qj) = 1. ■ 

It is not hard to see that a result similar to Lemma 5.1 can be established for con- 
jugate states tpij, which is stated below. Note that we can use the same C(z) and its 
multiplicative inverse B(z) in computations involving moments of A in the state ip as 
well as in its conjugate states ipij. 

Lemma 6.2. An X-valued power series lZ A (z) = 2^Li c n zn ~ x converging in the norm 
topology in a neighborhood of zero satisfies the equation (6.10) if and only if 

m 

(6-11) Yi Ti ViA h ni A ij h n 2 ■ ■ ■ &n fc _iAA fc ) = 

fc=l rii + ...+rife=m— k 
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for all m 2, where we assume that n±, . . . ,nk are non-negative integers and where 
= ^n=obnZ n+1 is the multiplicative inverse of C(z) = 1/z + lZ(z). 

Proof. The proof is similar to that of Lemma 5.1. ■ 

Theorem 6.1. Under the assumptions of Theorem 5.1, the series of the form 

(6.12) K A {z)=YjK i j(z) 

is an operatorial R-transform of the ipij- distribution of A it j for any (i,j) e J. 

Proof. Lemma 6.2 enables us to carry out the same approximation procedure for 
arbitrary arrays of strongly matricially free random variables in terms of operators of 
Toeplitz type as in the proof of Theorem 5.1, which leads us to the desired conclusion. ■ 

It turns out that if we require an operatorial R-transform to satisfy (5.6) for the 
state ip and (6.11) for all states ipij, then it is unique, and therefore it must be the 
matricial R-transform. In particular, in the Fock-space framework, these are the states 
defined by Q and the array Note that there is no need to study conjugate states 

of orders higher than two. On the Fock-space level (or, on the level of Hilbert space 
representations), the reason is that the action of strongly matricially free operators 
depends only on the type of the first vector in any simple tensor which defines a state. 

Theorem 6.2. Let Tl t j(z) = T it j{z)l it j for any (i,j) e J, where each Tij(z) is a power 
series converging in some neighborhood of zero. Then there exists a unique series 

(6.13) K A (z)=Y J T i , j {z) 

which is an operatorial R-transform of the distribution of A in the state if and of the 
distribution of A^j in the conjugate state ipij for any (i,j) e J. 

Proof. Existence of an operatorial R-transform of this form is established by Theorem 
6.1. Conditions (5.6) and (6.11) uniquely determine ip(b m ) and <Pj(b m ) for all m's and 
j e {1,2} since they are simple recursions. In the case of ip, we showed in Remark 5.3 
how to solve it, but the recursions for ip 1 and <p 2 are treated in a similar manner. Now, 
using 

i,j i,j i,j 

where pojections (qij) are given by Remark 2.1, we uniquely determine 7i,i +72,2, 71,2 + 
72,2, 72,i+7I,i, and thus also 72,i+7i,2, which implies that Ba(z) is uniquely determined, 
and consequently, its multiplicative inverse Ca(z) and thus TZa(z) in the orthogonal 
form 

n A {z) =Y J Q i , j {z) 

where Qij(z) = Qi,j(z)qi ! j for any (i,j) e J, with Qi/s being X- valued power series 
derived in the proof of Proposition 5.2, and where we identify q^/s with elements of T 
defined as in Remark 2.1. This completes the proof of uniqueness. ■ 
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7. CUMULANTS 

In this section we return to the study of the moments of the strongly matricially free 
convolution in the state ip, which can be identified with the collection of moments 

(7.1) MA m ) : m ^ 0}. 

Our results on the matricial R-transform indicate that there should exist a natural 
combinatorial approach which would allow us to express these moments in terms of free 
cumulants rather than to introduce new cumulants. Having in mind Theorem 3.2, this 
would also give an alternative approach to the c-free convolution [4] and computations 
of its moments without using c-free cumulants. 

If /i is the distribution of a random variable with moments {M At (m) : m ^ 0}, then 

(7-2) M M (m)= J] r,(|7r 1 |)...r M (K|) 

7r = {7ri,...,7r r }eA/ r C m 

where the numbers (r M (n)) n> i are the free cumulants of n [19]. 

Let us apply free cumulants to the combinatorics of the strongly matricially free 
convolution of (fJ-ij)- Let F n (n) be the set of all mappings / e i3(7r) — > {1,2, ... ,n} 
called colorings of the blocks B(ix) of n, where n e N. Then the pair (n, f) plays the 
role of a colored partition with blocks denoted 

(7.3) B(ir, f) = {(tti, /), (tt 2 , /),..., (vr r , /)} 

to which we shall assign free cumulants in a suitable 'matricial' way. 

In the case of the free convolution, if we are given distributions fii,fJ>2 whose free 
cumulants of order k are given by (rj(fc))i^^ 2 , then we can express the moments of 
A* = A*i IS A*2 as 

(7.4) M,(m)= J] r(7r 1} /)...r(7r r ,/) 

(T,/) = {(Tl,/),-,(Tr,/)}eVC m (2) 

where 

(7.5) r(ir k ,f)=r j (\ir k \) whenever f(n k )=j, 

and where MC m {2) denotes the set of non-crossing partitions of the set {1,2, . . . , m} 
which are colored by the set {1,2}. 

We would like to use a similar formalism to describe the strongly matricially free 
convolution. In contrast to (7.5), where cumulants are assigned to blocks independently 
of other blocks, we shall now make them depend on the colors of their outer blocks. 
For that purpose, instead of a tuple of cumulants of order k, we shall consider an array 
of free cumulants (rij(k)) for each k, where rij(k) is the free cumulant of order k of 
the distribution where e J. Clearly, this is in agreement with the matricial 
formalism on which our approach is based. 

We assign cumulants to colored blocks of (n,f) e MC m (2) as follows. If the block 
(ir k , f) and all its outer blocks are colored by j, we assign to it the diagonal cumulant 

(7-6) r(n k ,f) = rjJ (\n k \). 

On the other hand, if the block (7r k , f) is colored by i and it has an outer block colored 
by another color, then we assign to (n k , f) the off-diagonal cumulant 

(7-7) r(n k ,f) = r i:j (\n k \) 
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2 

n 

r[n,f]=r hl r 2tl r [x.f]='"i,2 r i,2 r KM= r i,i r 2,i r i,2 »' M [C,fc]=r lil r| jl r lj2 

Figure 1. Partitioned colored cumulants 

where j is the color of the deepest outer block of Ttk which has that property (j # i). In 
the case when some /z^- does not appear in the considered array (fiij), we set Tij(k) = 
for any k (one can also formally set fiij = So, which makes all cumulants vanish). Note 
that we do not require this outer block to be the nearest outer block of -k^ as in [12-13], 
where the combinatorics is slightly simpler since it is based on Gaussian operators, but 
we do require it to be its nearest outer block which is differently colored. The examples 
of Fig.l (especially the last one) show this feature in more detail. 

Therefore, if 7r fc is a block colored by j and all its outer blocks are colored by j, 
we assign to it the pair (J,j). If itu is colored by i and it has outer blocks wich are 
differently colored, then we assign to it the pair (i, j), where j is the color of the deepest 
blocks among those. We will then say that TCk is labelled by or respectively. 
The labelling induced by the coloring will be called J -admissible if e J. 

In order to eliminate labellings which are not J-admissible from our formulas, we will 
denote by AfC m (2, J) the subset of A/Yj m (2) consisting of those partitions which induce 
J-admissible labellings. For instance, if (1,2) ^ J, then all colored partitions in Fig.l 
but the first one are not J-admissible. 

Definition 7.1. Let /i be the distribution of A given by (3.1) and let (faj) be the 
corresponding array of distributions. By the partitioned colored cumulant corresponding 
to the colored partition (n, f) e MC m {2) we understand the product 

(7.8) r^,f]=r(n 1 J)...r(n r J) 

where tt = {tt\, 7r 2 , . . . , 7r r } and free cumulants are assigned to its colored blocks accord- 
ing to (7.6)-(7.7). 

As we already remarked in Section 4, it suffices to consider Toeplitz operators to com- 
pute the moments of EEkjA^j- We use these operators to derive the moment-cumulant 
formula given below. In contrast to the free or the c-free case, this formula is not the 
definition of cumulants since it expresses the moments of the sum of random variables 
in terms of free cumulants associated with marginal laws. 

Lemma 7.1. If fi is the ip -distribution of the sum A = X!jj a «,j of strongly matricially 
free random variables and fiij is the tpij -distribution ofa it j, then 

(7-9) M,(m)= r,[n,f] 

(n,f)eMCm(2,J) 

where the number r^[7r, f] is the partitioned colored cumulant corresponding to (tt, /). 



n n 



i i 

n n 



n n 
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Proof. Without loss of generality, we can choose (ajj) to be a two-dimensional array 
of Toeplitz operators of the form (4.1), where 

m 

f-A'h) !>.,(;•)(/*/ ; . 

k=l 

for some positive integer to, where we adopt the convention that = lij for any 

e J- Clearly, we obtain 

(p(A m )= ^(Sji-'-^jJ 

((_h,ji),-,(i m ,jm))eA m 

where A m is the subset of J m , for which there exist operators 

he{e ik , k ,(ei Jk r- l ;p^l}, l^k^m, 

such that 

(7.10) <p(bib 2 ...b m ) #0. 

We claim that the set of tuples (bi, b 2 , . . . , b m ) which satisfy these conditions is in 1-1 
correspondence with AfC m (2, J). 

Namely, to each (ir, f ) e AfC m (2, J) we assign a unique tuple (bi, b 2 , . . . , b m ) according 
to the following rules: 

(a) if k corresponds to the first leg of a p-block labelled by we assign to it 

h {it,)" '• 

(b) if k corresponds to any but the first leg of a p-block labelled by we assign 
to it b k = e id , 

where by a p-block we understand a block consisting of p elements. We claim that in 
this case (7.10) holds. In fact, if {k, . . . , k +p — 1} is a block for some k and p, we have 

(7.H) (^A+i,--- A+p-i) = ',,) 

where 

(7.12) = (i k ,jk) = ■■■ = (ik+p-i, jk+p-i)- 

Using (2.8), we can see that the associated product produces a power of a^j, which 
reduces the computation of the given moment to a moment of order s < m and we can 
apply a similar procedure to this reduced moment. The whole procedure finally gives 
a product of o^/s and thus (7.10) holds since a it j # for any e J by assumption. 

Conversely, suppose that a product bib 2 . . . b m of the considered type, with a nonvan- 
ishing moment in the state ip, is given. We would like to find a unique (rr, f) e AfC m (2, J) 
associated to this product according to the above rules. This will be done by induction 
with respect to to. First, observe that if to = 1, then the non- vanishing moment must 
be of type ip(ljj), where g J by Proposition 2.1(2), which corresponds to the 

partition consisting of a single 1-block colored by j with a J-admissible labelling 
If to > 1 and bk = ^i k> j k for all 1 ^ k ^ to, then we must have ik = jk and (jk,jk) G J 
for all k, which forces (n, f) to consist of 1-blocks only (they all have J-admissible 
labellings) and thus they must be colored by ii, i 2 , ■ ■ ■ , i m , respectively. 

If to > 1 and not all operators are units, there exists 1 ^ k ^ to such that (7.11) and 
(7.12) hold for some p ^ 1 since otherwise either all Vs would be creation operators, or 
all annihilation operators would be followed by other annihilation operators or 'wrong' 



26 



R. LENCZEWSKI 



creation operators, which would make the moment vanish. From among tuples of 
operators of this type, it is convenient to choose the one with largest k, which implies 
that bk+p-i = b m or bk +p is a creation operator with ik +p = j- Again, using the relation 
(2.8), we can see that such a product contributes a power of a it j. This reduces the 
moment to a non- vanishing moment of order s < m, to which there corresponds a unique 
7r' e MC s {2 ) J) with a J- admissible labelling by the inductive assumption. Now, there 
exists a unique ty e AfC m (2, J) obtained from it' by inserting the block {k, . . . , k + p — 1} 
colored by i in between k — 1 and k + p. Of course, if k + p — 1 = m, this block is 
separated from it' and we must have i = j since off-diagonal creation operators kill the 
vacuum, which uniquely determines the J-admissible labelling of the considered block. 
Otherwise, the inserted block is inner with respect to a block containing k + p and its 
J-admissible labelling is uniquely determined by the color of that block. 

Recall that the contribution from the product of operators given in (7.11), to which 
we associate a p-block, should be multiplied by the coefficient Sjj(p), which gives the 
contribution 

s m(p)«2 p-1) = r \j(p) 

from this block to the moment f(cii 1 ,j 1 cii 2 j 2 . . . a>i m ,j m )- The product of these contribu- 
tions gives the product of the cumulants corresponding to all blocks of (tt, f). Therefore, 
the contribution from the product of operators assigned to any (n, f) e MC m (2, J) is 

r(n l7 f) . ..r(ir r ,f), 

which completes the proof of the combinatorial formula for the moments of A. ■ 

Example 7.1. To colored non-crossing pair-partitions shown in Fig.l we assign par- 
titioned colored cumulants. For instance, we have r[n, /] = rfa, f)r(n 2 , f), where 
7Ti = {1,4} and tt 2 = {2,3}. Since all free cumulants involved are of order 2, we write 
rij instead of r i ^{2) ) which simplifies the notation. If we collect all colorings for the 
considered partitions, we obtain the sums of partitioned colored cumulants over all 
colorings: 

r M = r l,l( r l,l + r 2,l) + r 2,20"2,2 + ri )2 ), 

r [x] = r i,i( r i,i + r 2) if + r 2>2 (r 2>2 + r 1)2 f 
r K] = ^i,i( r i,i +r 2,i) 2 + r i,i r 2,i r ?,2 

2 / \2 2 

+ r 2,2l r 2,2 + r l,2) + r 2,2 r l,2 r 2,l 

Since £ is a partition of depth three, we can observe the influence of the definition 
of strongly materially free r.v. Namely, if blocks Ci = {1>8} an d C2 = {2,7} are 
differently colored, then blocks which are inner with respect to £2, namely (3 = {3,4} 
and (4 = {5,6} have to be colored differently than ( 2 due to the fact that we consider 
the case of strongly matricially free random variables. If we set r± t2 = r^i = r x and 
?"2,i = ^2,2 = f 2 (row-identically distributed square array), we obtain 

r[ir] = r\+ 2r 1 r 2 + r\, 

r [x] = r i + 3r^r 2 + 3rir\ + r 2 , 

r K] = r\ + Srfr 2 + 2r\r\ + 3nr 3 2 + r\, 

which is the contribution from partitions n, x an d C t° the moments of [i\ El A*2- m 
turn, if we set r 2) i = r 2j 2 = r 2 , r^i = r\ and r 12 = (row-identically distributed 
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lower-triangular array), we obtain 

r[n] = r\ + r\r 2 + r 2 , 

r[x] = ?"i + 7"ir 2 + r x r\ + r|, 

r[C] = + 2r\r 2 + r\rl + 4, 
which gives the contribution from ir, x an d C to the moments of ji\ o ^ 2 . 

Example 7.2. The lowest order moments of A are expressed in terms of free cumulants 
of the distributions fi^j as follows: 

M„(l) = r 1)1 (l)+r 2 , 2 (l), 

M M (2) = r lil (2)+r 2 , 2 (2) + (r 1 , 1 (l)+r 2 , 2 (l)) 2 , 

M M (3) = n,i(3) + r 2l2 (3) + 2(n,i(2) + r 2 , 2 (2))(r 1 , 1 (l) + r 2;2 (l)) 

+ r 1;1 (2)(r 1)1 (l) + r 2jl (l)) + r 2)2 (2)(r 2>2 (l) + r 1)2 (l)) 

+ (r 1;1 (l)+r 2 , 2 (l)) 3 . 

As in the case of moments of arbitrary orders, these moments agree with the moments 
of fj,i El A*2 if the array (ajj) is square and row-identically distributed. In this case 
the moments of A can be expressed in terms of free cumulants of fii EB A*2- m turn, 
if that array is lower-triangular and row-identically distributed, these moments agree 
with those of o // 2 . However, as already M M (3) demonstrates, the partitioned 
colored cumulants which appear in our moment-cumulant formula cannot be expressed 
in terms of free cumulants of free convolutions of distributions. Roughly speaking, this 
effect shows how much the shape of the lower-triangular array affects the additivity of 
the considered transforms. 

8. Cauchy transforms 

Continuing the considerations of Section 7, we shall study the Cauchy transforms of 
the 'commutative' (^-distribution of A and express them in terms of their R-transforms. 
One should remember that these are Cauchy transforms with complex arguments rather 
than those with operatorial arguments which provided the right framework for the 
linearization property of the R-transform. 

As a consequence of Proposition 4.1 and Lemma 7.1, the (/^-distribution of A is 
uniquely determined by the array of free cumulants of distributions (^j). The way 
we defined partitioned colored cumulants indicates that the moments of A under <p 
exhibit the property of 'diagonal subordination', which is reflected by the fact that all 
covering blocks of non-crossing partitions which appear in the formula for the moments 
of A in terms of free cumulants are labelled by diagonal pairs where j e {1,2}, 

and all inner blocks have colorings which are 'matricially subordinate' to their closest 
outer blocks. 

More precisely, the first level (counting from the top) of each non-crossing partition 
(i.e. all covering blocks, including singletons) is labelled by some but all blocks 

which are inner with respect to a given block iTk labelled by and which form 

nearest inner-outer pairs with that block (thus, all blocks at the second level counting 
from the top) are 'subordinate' to 7tk and the corresponding coloring in the sense that 
its labelling must be for some i. A similar 'subordination' holds on the remaining 
levels. 
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Figure 2. Binary subordination of measures 



Therefore, it is natural to introduce a mapping which assigns to each distribution 
Hij the variable A it j associated with the array (fJ>*j) which is 'subordinate' to fj, it j (to 
the top covering block we assign r^j). 

Definition 8.1. Let (//*,_,-) be an array of distributions. By the subordinate arrays of 
distributions we understand arrays of the form 

(8.i) a 1= ( w ^ y ti 2 =( ^ ^ y 



(8.2) £ 1= ( ^ ^ y £ 2 =( ^ ^ ^ . 

' V ^2,1 A*2,l / ' V ^2,1 A*2,2 7 

The corresponding Cauchy transform will be called subordinate Cauchy transforms and 
will be denoted by G* -, respectively. 

Figure 2 shows how to interpret the subordination of distributions exhibited by these 
equations. It has the form of a binary tree, where the part of the tree that lies below 
any /i^ shows which measures have to be replaced in the original array in order to 
produce a formula for the corresponding moments of strongly matricially free random 
variables. 

More explicitly, these are Cauchy transforms of the sums of strongly matricially free 
random variables with distributions given by (/-f*^)- It is easy to see that they satisfy 
the recursion given by equations 



with 1=2 and 2 = 1, where K itj (z) = R i j(G* j (z)) for any i,j e {1,2}. 
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Theorem 8.1. Under the assumptions of Lemma 7.1, the Cauchy transform of the 
distribution of A in the state ip is given by 

z -Zij R 3A G jA z )) 

where (G*j) is the associated array of subordinate Cauchy transforms. 

Proof. The general combinatorics of the proof reminds that for conditional freeness. 
The most important difference is that we use free cumulants instead of c-free cumulants. 
Let 

M A (z) = z-'GAiz- 1 ) 

be the moment generating function associated with the (^-distribution of A. Then, by 
Lemma 7.1, we have 

oo 

M A (z)-l = Y,z n 2 r[n,f]. 

n=l (n,f)eATC c n 

On the other hand, if our assertion is supposed to hold, we must have 

2 

M A (z)-l = J]K jd (z- 1 )zM A (z) 

3 = 1 

where 

oo oo 

Kjjiz- 1 ) = S^'H 2 {zMl j {z)r...{zMl j {z)T^ 

m=l ni,...,n m _i=0 

which leads to the equation 

2 oo oo 

m a {z) - 1 = 2 2 r «H 2 m *>o • • • ^("-iww^ 1 ^- 

j = lm=l ni,...,n m =0 

Therefore, the coefficient standing by z n in this expression is given by 

2 n 

C A(n) = 22 2 r j,j( m ) m j,j( n l) ■ ■ ■ m *jK-l) m -4( n ™) 

j = l m=l m + . . .+n m =n—m 

We need to compare this number with m A (n) expressed in terms of free cumulants: 

m A (n) = 2 r [ n J] 

If we express each moment m*j(rik) in the above formula for c A (n) in terms of free 
cumulants, we can argue that we obtain the same expression as that for m A (n). In fact, 
it suffices to interpret the product 

r j,j( m ) m j,j( n i) ■ ■ ■ mlj(n m -i)m A (n m ) 

as follows. The cumulant r^j{m) is assigned to this block of some (it, f) e MC c n which 
contains 1. This block is assumed to have m legs and to all of these legs but the first 
one we assign 'subordinate' moments m*j(ni), . . . , m*j(n m _i) (each of these moments 
is a sum of paritioned colored cumulants associated with the partitions of the subin- 
terval lying between the given leg and its left neighbor). Here, numbers n 1; . . . , n m _i 
are cardinalities of these subintervals. This part of the product corresponds to the first 
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covering block and its inner blocks. It remains to assign m J 4(n m ) to the remaining 
covering blocks and its inner blocks which has the same general form as m^ra) except 
that the cardinality of the considered interval is smaller. It can be seen that each par- 
titioned cumulant is obtained in this fashion and it is obtained precisely once since the 
numbers j, m, cardinalities ni, . . . , n m _i of the above-mentioned subintervals, together 
with m — 1 colored non-crossing partitions of these subintervals, and a fixed product 
of free cumulants taken from the cumulant expression for m,A{n m ) determine a unique 
partitioned cumulant of (n, f). This completes the proof. ■ 



Corollary 8.1. Under the assumptions of Theorem 8.1, the Cauchy transforms of the 
following convolutions are expressed in terms of R-transforms by formulas given below: 

(1) free convolution 

(2) monotone convolution 

1 

G>ll>|li2 ( z ) = 



(3) boolean convolution 

1 



™ V ~ y z-R^G^-R^G^z)) 

(4) s-free convolution 

(5) orthogonal convolution 



/i 



h3 



z R^ (G^ ll>jU2 (z)) 

Proof. We shall present the proof of the first two formulas. In the case of free 
convolution, we have 

since distributions are row-identical and therefore the array //* • agrees with the original 
array (faj). Hence, G* ■ = for any and thus (8.3) easily gives the formula for 
the free convolution. In turn, in the case of monotone convolution, the array is lower- 
triangular and ri^(k) = for any k. This implies that all blocks lying under blocks 
colored by 2 must also be colored by 2. Therefore, /i| 2 reduces to the one-dimensional 
array with ^2,2 = A*2, whereas 

since blocks lying under the covering blocks colored by 1 can be colored by 1 or 2. 
However, if such a block is colored by 2, then the deepest block which is differently 
colored and is outer with respect to that block, must be colored by 1. This gives 
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= C Mll>M2 and G 22 = G^, which gives the formula for the monotone convolution. 
The proofs of the remaining formulas are similar and are left to the reader. ■ 

Remark 8.1. The equations of Corollary 8.1 can also be proved without using Theorem 
8.1. For instance, using formulas 

i^ioiAj = F Vl o F U2 and F v {z) = z - R u (G u (z)), 

where F u is the reciprocal Cauchy transform of v, we obtain 

F^ 2 (z) = F^z) - R^G^F^z)) 

= z — R fl2 (G f , 2 (z)) — Rn 1 (Gf 11 c IJl2 (z)), 

which is equivalent to equation (2). In a similar way, one can prove the remaining 
equations. 

Nevertheless, Theorem 8.1 gives a formula which relates Cauchy transforms of various 
convolutions to R-transforms of \i\ and \x 2 in a unified manner. Let us also remark that 
this formula can also be viewed as an approximation for Cauchy transforms of various 
convolutions. In that respect, it is similar to formulas for the Cauchy transform of the 
free convolution of compactly supported measures given in [10]. 

Example 8.1. Let us consider the example of a two-by-two square array in which 
the diagonal distributions are semicircle laws whereas the off-diagonal ones are point 
masses. Thus, let R\^{z) = az, #2,2(2) = dz and R lj2 (z) = b, R 2jl (z) = c. Then, using 
(8.3), we obtain 

G " {Z) = z - aG* 1:1 (z) - dG* 2:2 (zy 
where fi = EEkjAtjj and the diagonal subordinate Cauchy transforms correspond to 
shifted semicircle laws since they satisfy the equations 

Gii( z ) = — -. z and G 22 (z) = : — — . . . 

1,1 v ; z-c-aGl tl (z) 2 ' 2V 1 z-b-dG\ 2 {z) 

In the general case, the analytic formula for G^ is rather complicated, but if a = d and 
c = b, then G* 1 = G 2 2 and thus 

1 6-V(&-^) 2 -4o 
A } z-2aGl 1 (z) 4a + 2bz-z 2 ' 

and therefore /i is the free Meixner distribution with the continuous density 



7r(4a + 2bx — x 



2^ 



on the interval [b — b + 2*J~a\ and two atoms at b + \/b 2 + 4a. 

Note that the computations of the Cauchy transforms of strongly matricially free 
convolutions differ from those in [12,13], where certain examples of matricially free (but 
not strongly matricially free) convolutions arising in the central limit theorem were 
studied. 
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